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Abstract 

The effect of retaining the material dispersion terms in the nonlinear coupled mode equations (NLCME) 
that describe light propagation in fiber Bragg gratings is analyzed. It is found that dispersion is responsible 
for new instabilities of the uniform states and gives rise to new complex spatio temporal dynamics that is 
not captured by the standard NLCME formulation. A detailed analysis of the effect of dispersion on the 
linear stability characteristics of the uniform solutions is presented and some numerical integrations of the 
NLCME with dispersion are also performed in order to corroborate the theoretical results. 
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I. INTRODUCTION 



The nonlinear coupled mode equations (NLCME) are the envelope equations currently 
used to study the weakly nonlinear dynamics of light propagation in fiber Bragg gratings. 
These are universal equations in the sense that they describe the large scale pattern forma- 
tion of any conservative, spatially extended propagative system with a weak periodic spatial 
structure. Material dispersion effects are systematically neglected in the derivation that leads 
to the NLCME. This paper presents a study of the effect of retaining the dispersive terms 
on the stability of the uniform states and it concludes that these terms should be considered 
because they give rise to new instabilities and dynamics that are not captured by the NL- 
CME. The NLCME fail, in general, to predict the system evolution and the dispersive terms 
produce the onset of short scale perturbations that evolve and spread over the domain in 
the NLCME time scale. Numerical results of the integration of the NLCME with dispersive 
terms are also provided to corroborate the theoretical linear stability results. 

Fiber Bragg gratings (FBG), i.e., optical fibers with a periodic variation of the refractive index 
along its length, exhibit frequency gaps where light resonates with the periodic structure of the 
fiber and propagation is not permitted. This frequency dependent reflection characteristic together 
with the recent developments in grating fabrication techniques has allowed for successful applica- 
tion of FBG in optical communications and fiber sensing, see e.g. jl] for a recent review. Apart 
from this purely linear effect, the combination of the nonlinearity of the fiber that comes into play 
at sufficiently high intensities with the spatially distributed reflection produced by the grating gives 
rise to localized structures known as gap solitons (see |2, 3, ^ and references therein). These gap 
solitons are solitary waves that, at least in theory, can propagate through the FBG at any speed 
between zero and the speed of light of the fiber without grating. This possibility of capturing 
light (gap soliton propagation at zero speed) is currently a topic of great interest because of its 
promising future applications in the production of devices for optical storage of information. 
Light propagation in a nonlinear FBG is usually analyzed using the so-called nonlinear coupled 



mode equations (NLCME) 



A+ = VgA+ +iAnA- +igA+{a\A+\^ + \A-\^), (1) 
At = -VgA- + iAnA+ + igA-{a\A-\^ + \A+\^), (2) 

which prescribe the time (t) evolution along the fiber (x) of the complex envelopes ^4+ and A~ of 
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the wavetrains that resonate with the grating. All variables are assumed to have been previously 
nondimensionalized in a similar way as in [|61], that is, space has been made nondimensional with 
the grating period Ag, time with Ag/c (c is the vacuum speed of light), the grating amplitude with 
the mean value of the refractive index, and the fields with the characteristic size that results from 
the Kerr nonlinearity. The actual electric field inside the FBG can be written, in first approxima- 
tion, as the superposition of two modulated, counterpropagating wavetrains 

E ~ A+{x, t) + A-{x, t)e-^™-+^"«* + ex. + ■ ■ • , (3) 

with wavenumber vr and frequency ujq = uj{k = vr), which is given by the dispersion relation 
of the fiber without grating u{k) (c.c. stands for the complex conjugate). The wavelength of the 
wavetrains, Aq = 2, is twice the period of the grating, that is assumed to have a spatial profile of 
the form cos(27rx). This situation can be seen as a spatial parametric resonance that couples both 
wavetrains: the interaction of one of the wavetrains with the grating gives rise to a resonant term 
that excites the other wavetrain, 

cos(27rx)A+e'™'+^^«* ~ (e'^"^ + e-i2-^)/l+ei™-+i^ot ^ ^+g-i-x+ia;ot j-gsonant), 
and viceversa. 

The NLCME take into account the effect of (i) transport at the group velocity, which is given by 
Vg = duj / dk dXk = TT, (ii) coupling of the wavetrains through the spatially extended reflection pro- 
duced by the grating, with An proportional to the strength of the grating, and (iii) cubic nonlinear 
interaction, where a = | as a result of the cubic Kerr nonlinearity of the original physical problem 
for the electric field 121]. The NLCME are also used to describe the evolution of Bose-Einstein 
condensates in optical lattices 13, [SD and, in general, the NLCME can be regarded as a normal 
form that applies to any dissipationless propagative system, extended in one spatial dimension, 
reflection and translation invariant, and with a spatial parametric periodic variation. 

The derivation of the NLCME is carried out using a standard multiple scales procedure that 
requires to assume small amplitudes in ([3]) with slow spatial and temporal dependence (slow as 
compared with the wavelength and temporal period of the basic wavetrains), i.e., 

•■• < < < lA^I < 1, ■ ■■ < < < 1, (4) 
and small grating intensity. 

An < L (5) 
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(7) 



The rest of the parameters: ujq, Vg and g, are all considered order one quantities. A detailed 
derivation of the NLCME from the ID Maxwell-Lorentz equations can be found in [6]. In the pro- 
cess of derivation of the NLCME material dispersion terms (~ i^^^^) are systematically neglected 
against transport terms (~ A^). Notice that this is not the case in the derivation of the nonlinear 
Schrodinger equation (NLS) for a single modulated wavetrain propagating in an uniform fiber, 
where both effects are considered because the transport effect can be easily accounted for just by 
using a moving coordinate frame that travels with the group velocity If we retain the material 
dispersion terms, then the resulting equations for A^ read 

Af = VgA+ + idA+^ + i/\nA~ + igA 
At = -VgA- + idA-^ + iAnA+ + igA 

where d = —d'^uj/dk'^ at /c = vr is the dispersion coefficient that is, generically, an order one 
quantity. Note that the coefficients of the linear terms of the above envelope equations can be 
easily obtained from the coefficients of the power expansion of the dispersion relation uj{h:) at 
k = ±TT (see e.g. the review lioll '). 

In the above equations the transport terms cannot be removed using a moving coordinate trans- 
formation and the two effects, transport and dispersion, with different asymptotic orders, see eq. 
must be considered simultaneously and give rise to two different slow spatial scales: 

a) A transport scale, (5xtrans ~ 1/ An ^ 1, which comes from the balance between the transport 
along the fiber at the group velocity and the continuous spatial reflection produced by the grating 
and is the slow spatial scale used in the NLCME derivation flflH]. 

b) And a dispersive scale, 5xdisp ~ 1/ VAn ^ 1, which results from the balance between disper- 
sion and grating reflection. This dispersive scale is small as compared with the transport scale 
but still large as compared with the basic wavetrain wavelength, i.e., 

and corresponds to the standard dispersive scaling used in the derivation of the NLS equation 
10. 

If we look for states exhibiting only transport scales then the dispersion terms in equations ©-(IT]) 
can be neglected in first approximation and these transport dominated solutions are thus well de- 
scribed by the NLCME. But, in general, the question of whether the system will develop dispersive 
scales or not is a stability question and its answer is not known a priori. 



The purpose of the present paper is to demonstrate that the material dispersion terms in eqs. 
©-(IT]) should not, in general, be neglected because the dispersive scales can be unstable and grow 
and spread all over the domain. In other words, we will show that stable solutions of the NLCME 
can be unstable in the context of the more general problem given by eqs. ©-(jT]), and can develop 
dispersive scales that give rise to very complicated spatio-temporal states that are not included in 
the NLCME description. 

To this end, we will consider equations ©-dV]) with simple periodic boundary conditions, 

A'^ix + L,t) = A^{x,t), withL>l, 

which correspond to a ring shaped FBG with length L containing an even number of grating 
periods (like the fiber-loops considered in the simulations of multiple gap soliton collisions in 



llllh ). and we will focus on the effect of dispersion on the stability of the simplest solutions of the 



NLCME, namely, the solutions with constant uniform modulus. Equations ©-(IT]) can be further 



simplified if we rescale them with the transport scales commonly used for the NLCME 



iQE, 



t = tVg/L, x = x/L, = \J gL/vgA , and K = AnL/vg, (8) 

where the scaled grating strength k ~ 1, i.e., the FBG length is of the order of the transport scale 
^a;trans- Notc that g can always be made positive by making, if necessary, the change of variables 
A^ -^A^ and taking complex conjugates in eqs. ©-dT]). The resulting nonlinear coupled mode 
equations with dispersion (NLCMEd), after dropping tildes, take the form 

At = A+ + i£A+ +ifi:A-+iA+(a|A+|2 + |A-|2), (9) 
Ai = -A- + ieA-^ + iKA+ + iA-{a\A-\^ + \A+\^), (10) 
A^{x + l,t) = A^{x,t). (11) 

where the scaled dispersion e = d/{Lvg) <^ 1 and k ~ 1, as mentioned above. With this new 
scaling the transport and dispersive scales are now (5xtrans ~ 1 and 5xdisp ~ ^ 1' ^^'^ the 

wavelength of the basic wavetrains in dS]) is ~ |£:| ^ 1. The small dispersion term e <^ 1 in 
the above equations is essential because it takes into account the fact that we are simultaneously 
considering two effects, transport and dispersion, that have different asymptotic orders due to the 
original slow modulation assumption (H]). There are some previous works OHQ that study 
the structure and stability of the localized solutions of the NLCMEd (I9l)-([T0]). but they all consider 
£ ~ 1 and do not study the small dispersion limit e <^ 1, which is the only physically relevant 



regime from the point of view of the slow envelope description of the weakly nonlinear light 
propagation in FBG. 

The remainder of this paper is organized as follows. The family of uniform solutions of the 
NLCME is briefly described in the next section and its linear stability characteristics against per- 
turbations with wavelength A ~ 1 (transport scale) and A ~ \/W\ ^ 1 (dispersive scale) are 
obtained in sections 3 and 4, respectively. These results indicate that the dispersive scales can 
destabilize the uniform solutions for parameter values where the NLCME predict stability. In the 
final section, numerical simulations of the NLCMEd with small dispersion are used to corroborate 
the theoretical stability predictions and some concluding remarks are also drawn. 



II. CONTINUOUS WAVES 



The constant modulus solutions of the NLCME, also known as continuous waves (CW), can be 



written as y. 



= pcos^e'"*+^™^, 
= psin0e'"*+''"^, 

where p > is the power flowing through the grating, 9 e] — |,0[U]0, |[ measures the ratio 
between the two counterpropagating wavetrains, and a and m are given by 

K 1 + CT 

K 1 — a 2\ 
"^ = (^ + ^^^^°^'^- 

According to eq. ([3]), the associated electric field inside the grating consists of the superposi- 
tion of two uniform wavetrains, and a and m represent small corrections to the frequency and 
wavenumber of the wavetrains. The resulting pattern resembles a standing wave for 9 near ±^ 
and a traveling wave for 9 close to and ±|. 

In the linear limit of small light intensity, p ^ 0, the relation between the frequency, a, and the 
wavenumber, m, of the CW takes the form 

^2 

= ^), 

which is plotted in Fig. 1 and shows a frequency gap of size 2k, where light propagation is not 
allowed (this is a typical resonance effect in wave propagation in periodic media). As the power 




Figure I: a — m plot of the CW in the linear regime p — > 0. 



b) 




Figure 2: a — m plot of the CW for a) p < pc and h) p > pc (a = 

is increased, the nonlinearity of the problem comes into play and the relation between a and m 
becomes more involved 
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There are still two CW for any given m if < = 2k/{1 — a), see Fig. 2a, while for p^ > p^ a 
region of higher multiplicity of solutions (up to 4) develops near m = I4], as it can be appreciated 
from Fig. 2b . Notice that the frequency gap persists for p > and it remains of size 2k but it is 
shifted upwards by the nonlinear effects. 

In the rest of this paper the family of the CW is going to be represented in the (9, p^) plane 
shown in Fig. 3. The CW inside the region marked by thick line 

= = -71 T- (oa\ (12) 

(1 — a) sm(2t^) 



K 




e 



TT TT TT TT 

2 ~4 4 2 



Figure 3: {9, p^) representation of the CW {a = \). 

are those along the loop in Fig. 2b. The wavenumber m is positive (negative) in the white (shaded) 
areas and the thin lines correspond to the values of the wavenumber compatible with the periodicity 
boundary conditions (fTTI) . m = 27m with n E Z (arrows indicate the m increasing direction). Note 
that the 6* = ±1 symmetry is due to the fact that the CW associated with 9 > (9 < 0) and j - 9 
(— I — 6*) are essentially the same after applying the spatial reflection symmetry of the system 
X —X and A'^ ^ A^, 

III. CW STABILITY: TRANSPORT SCALE PERTURBATIONS 

The CW are approximate solutions (up to order e corrections) of the NLCMEd and its infinites- 
imal perturbations 

= A+„(l + a+), A~ = A^^{l + a~), with |a^| < 1, 

evolve according to the linearized version of the NLCMEd 

— a'^ = iK{a~ — a"*") tan 9 + iap^ cos^ 9(a'^ + a+) + ip^ sin^ 9{a~ + a~) + iea^^, (13) 
+ a~ = iK(a^ — a~)/ tan 9 + iap^ sin^ 9{a~ + a") + ip^ cos^ 9{a~^ + a+) + iea~^, (14) 
a^{x + l,t) = a^{x,t), (15) 
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which, by means of the spatial Fourier expansion 



k=—oo 

can be turned into a system of ordinary differential equations of the form 

da 



= i(27r/c)a^ + iK(a^ - a^) tan 6^ + iap^ cos^ 6^(0^+ a+^) 



+ipsin^ e{ai. + a^,^) - ie{2TrkyaJ, (16) 

dal 



''k 



-i{27ik)aj^ + iK{a'^ — a/^)/ tanO + lap sin 6'(a^ + a_ 



dt 

+ip^ cos^ 6^(0^ + atk) - ie{27rkfa^, (17) 

where k E 7j indicates the wavenumber of the perturbation. 

Perturbations with wavelength of the order of the transport scale have wavenumbers A; ~ 1 and 
for them the dispersion terms ie{27rk)'^a'^ in the equations above are small and can be neglected in 



first approximation. The resulting equations for together with those corresponding to aZf, form 
a linear quartet uncoupled from the rest whose solutions are of the form 



{4{t),at,{t)) = {al^,a'iy-\ (18) 



{lo^ - {2'Kkff - 2k\uj'' - {2'Kkf) + Akp \ ((27rfc)^(l + a) - uj\l - a)) 



where uj is given by the following fourth order polynomial with real coefficients 

tan^ 
1 + tan^ e 

-K^ tan^ eiuj + {2Txk)f ^(w - (27rA;))2 = 0. (19) 

tan B 

This polynomial results from the NLCME, that is, when dispersion is not considered, and it was 
previously analyzed in lilSll for some particular values of (? and Q. In this section we will complete 



these results and give a detailed description of the CW instability regions given by (fT9l) . 

For the particular case of uniform perturbations. A; = 0, the solutions of ( fT9l) can be calculated 
explicitly 

= (double) and 



sii?(2^) 



^ = ±a/ ■ + 2Kp2(l - o) sin(2^), 



and we have instability, i.e., uj with negative imaginary part, when 




Figure 4: CW stability properties (cr = Shading indicates instability and the thin lines correspond to 
ki = 0.25 + 0.25n, with n = 0,1,2, ... . The CW above the ki = 2tt/k line are unstable, see eq. (|23] ). 

This instability region lies inside the higher multiplicity region given by eq. (fT2l) . see Fig. 4, and 
it destabilizes the CW on the upper part of the loop in Fig. 2b. 

In the limit of large wavenumbers, k ±oo, the solutions of eq. (fT9l ) can be expanded as 

u = ujq + LJi + . . . , with IljqI ^ l^il ^ • • • ■ 

The first order term cuq ~ A: ^ 1 satisfies the equation 

{col - i^^kff = 0, 

whose solutions are of the form uJo= ± 27rk and represent pure oscillations. The next order correc- 
tion cji ~ 1 is given by 

Q o If) o ijSiii 

LOi = hi tan — zanp 



1 + tan^ 9 

and the CW are therefore unstable inside the region defined by the conditions 

9 9 K tan^^ 6 

P >Poo= ■ .^ny (21) 

which is plotted in Fig. 4. Notice that, for the CW in this region, all perturbations with wavenum- 
ber above a given one are unstable. This unphysical behavior (amplification of perturbations far 
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away from the Bragg resonance) was already detected in 11511 where it was suggested that it could 
be fixed by the addition dispersion. This affirmation is true, as it will be seen in the next section 
where the stability of perturbations with high wavenumbers is analyzed. 

The localization of the unstable perturbations with finite nonzero wavenumber < |A;| < oo 
requires to numerically explore the roots of (fT9l) and can be simplified by taking into account the 
invariance of eq. (fT9l) under the transformations 



{uj, k) ^ (—to, —k) and (co',tan6') — (— cj, l/tan6') 

that allows us to restrict the search to the parameter range 

k>0 and ^€[--,-1. 

4 4 

Furthermore, if we use the scaled variables 

u = u/k, k = {2Trk)/K and = / k, (22) 
then the strength of the grating, k > 0, can be absorbed and eq. (fT9l) can be rewritten as 

(^2 _ ^2)2 _ 2{Co' - P) + 4p2-i^^(P(l + a)-u\l- a)) 
- tan2 eiCj + ~kY X^iCj - kf = 0. 

The onset of complex roots (instability) takes place when both the equation above and its derivative 
with respect to u vanish. From these two conditions the values of 9 and fp can be numerically 
computed for any given u and k, and the following results are obtained: 

• If > then there are no more instabilities apart from the one already obtained for k ±oo 
dZB, see Fig. 4. 

• For each negative 9 and for each value of = /n there is a critical scaled wavenumber 
ki such that the perturbations with < A; < are unstable and those with A; > A;i are stable 
(constant ki curves are plotted as thin solid lines in Fig. 4). A CW will be then unstable if 
the perturbation mode with lowest wavenumber (k = 1) is unstable, that is, if 

— < ~ki, (23) 

K 

recall the definition of k in (|22l) . For a given value of k, all the CW above the line ki = 2tx / k 
in Fig. 4 are unstable. 
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• The CW near the = pi line exhibit a small interval of stable scaled wavenumbers k e 
[0, with ^2 < kl, but, as it has been checked numerically, this does not modify the 
instability region defined above. 

The compete stability characteristics of the CW against transport scale perturbations (which are 
accurately described by the NLCME) are summarized in Fig. 4. There are three instability bound- 
aries that mark the destabilization of the uniform perturbations (|20l) . high wavenumber perturba- 
tions (|2TI) and perturbations with wavenumber k = 1 (1231) . The two first instability conditions 
are not sensitive to the FBG length and depend only on the combination p'^/n while the third one 
depends also on k and it accounts for the stabilizing effect of the finite size of the domain; in an 
infinite FBG (k oo, see the last eq. in ([8])) all CW with 9 < are unstable due to (l23l) . 



IV. CW STABILITY: DISPERSIVE SCALE PERTURBATIONS 

The CW perturbations with wavelength of the order of the dispersive scale have large wavenum- 
bers, k ~ l/A/|e| ^ 1, and for them the dispersion terms in eqs. (fT6l)-(fT7l) can not be ne- 
glected. In order to study the stability of these perturbations, we define the scaled wavenumber 
K = {2Txk) a/I^I ~ 1 and expand the solution of eqs. (fT6l)-(fT7]) in the form 

= o-Koi^^ ^) + VWl^Kiit, T) + ■ ■ ■ , a]^ = aJ^Q{t, T) + ^/\e\aK^{t, T) + ■ ■ ■ , 
where T = t/ >/\e\ is a fast time scale. At first order we obtain 

^"^0 _ iKa+ - 
da 



+ iKa-^, = 0, 



whose general solution is given by 



Thus, in the fast time scale T, the group velocity term dominates and these narrow dispersive 
perturbations simply travel in opposite directions. At next order we obtain 

d(x^-\ . ^ r dA. 



^i-ii^a+i = [-^^-i(Ktane + ^K2)A+o + iap2cos2^(A+o + A+^J]e'^^ 



dT ' dt ' \e 



-[iKtanM^o + ip2 sin^ e{Aj^f^ + AZko)]^''^^ 



^^ + ^Ka],, = h^-i(«:/tan^ + -i^2)yl-o + iap2sin2e(A7,o + A:,,o)]e-'^^ 



+ [iK/tan^A+o + ip^ cos^ ^(A+q + A+^o)]e 
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iKT 



and the higher order correction (a^^, aj^-^) will then remain bounded in the fast time scale T if the 
following equations are satisfied 



dt 



i(Ktan0 + —K')A+Q + iap' cos' e{A+Q + A+^q)' (24) 



— ^ = -i(fi:/tan0 + ^/^')A^o + i^P' sin^ ^(A^^ + AI^q), (25) 
which give the evolution of A^q in the slow time scale t. The left (+) and right (~) propagating per- 



turbations are not coupled, and the equations above together with the corresponding ones for A_j^q 



and ^_^'o form two linear systems with constant coefficients whose solutions are proportional to 



= ±J{Ktan9 + -^K^){2ap'^cos'^9-{Ktan9 + -^K'^)) and (26) 



n- = ±^(/€/tan^ + ^A'2)(2o-p2sin2^-(K/tan^+ I^K2)). (27) 
The dispersive perturbations are then unstable if the following conditions are satisfied 

< K tan ^ + -^^K^ < 2a cos^ 9, 

^ -'2 ^ o^„2 „• 2 , 



< k/ tan ^ + — K' < 2crp2 sin' 9, 



\e\ 

that can be expressed as 

P >pi = --^(Ktan^ + T-K^) with tan^ > — — , (28) 

asm(26') \e\ \e\ k 

tan-ig e^, with tan-i^>^ — . (29) 

(Tsin(2^) l^l \e\ k 

Notice that if we set K"^ = in the equations above then the instability condition given by eq. (l2T]) 
is recovered, in other words, the limit k ±oo for the /c ~ 1 (transport scale) regime matches 
with the limit K — > for the /c ~ 1/ ^ 1 (dispersive scale) regime. 

The above instability conditions depend on the parameter > and define two families 
of regions in the {9, p^) plane, which are represented in Fig. 5 for e > and several values of 
K"^ > (note that expression ( |28] ) is identical to (|29l ) if we change tan 6^ by 1/ tan 9 and therefore 
the instability region defined by the condition is simply the symmetric of that defined by 

TV 

around the vertical axes 9 = A variation of the wavenumber from ktok + 1 results in a very 



small increment of the scaled wavenumber AK ~ ^ 1, so, in first approximation, we have 
to allow K'^ to vary continuously between and +oo in (l28])-(|29l). The resulting final dispersive 
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instability region is then given by the union of all these K'^ dependent regions. All CW with 9 
negative are thus rendered unstable for e > (see Fig. 5) and, in a completely similar way, if e < 
then all the CW with positive 6 are unstable (see Fig. 6). 




TV TV Q TV TV 

2 ~4 4 2 



Figure 5: Dispersive instability regions defined by eqs. (1281 ) and ( [29l ) for e > and three different values of 
K^, with arrows pointing towards increasing direction and shading indicating instability (a = 

Therefore, the CW stability diagram represented in Fig.4 has to be complemented with the 
following dispersive instability criterion: 




Vr TT TV TV 

2 ~4 4 2 



Figure 6: Dispersive instability regions defined by eqs. ( [28] ) and ( [291 ) for e < and three different values of 
K^, with arrows pointing towards increasing direction and shading indicating instability {a = |). 
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k - 1 



k ~ 1 



A'~l (fc ~ l/v/|£| > 1) 



A'~l (fc ~ l/v/|£| > 1) 



Figure 7: Sketch of the maximum growth rate of the transport and dispersive perturbations as a function of 
the wavenumber 



which states that all the CW along the lower (upper) branch in Fig. 2 are unstable if e is positive 
(negative). 



As it was advanced in [|l5n . the material dispersion corrects the unphysical high wavenumber 
instability found in the previous section, see eqs. (|26l) and (1271 ) and the upper plot in Fig.7. But 
it also gives rise to new instabilities (lower plot of Fig.7) that are not detected if the NLCME 
description is used and that are present independently of the dispersion sign: according to (l30l) . 
positive (negative) dispersion is responsible for the destabilization of NLCME stable CW with 
6* < (6* > 0) in the lower part of Fig. 2. These high wavenumber dispersive perturbations travel 
with the group velocity and their growth rate remain of order unity as £ — 0, see eqs. (l26l)-(IT7l). 
i.e., they evolve in the transport time scale t ~ 1. 



V. DISCUSSION AND CONCLUSIONS 

As a final check of the theoretical dispersive instability results obtained in the previous section, 
several numerical simulations of the NLCMEd have been performed for the two cases shown in 
Fig. 8 with different dispersion coefficients. The numerical method uses a Fourier series in space 
with A^^Fourier modcs and a 4th order Runge-Kutta scheme for the time integration of the resulting 
system of ODEs, with the linear diagonal terms integrated implicitly and the nonlinear terms 
computed in physical space using the 2/3 rule to remove the aliasing terms ll^ (typical required 
resolutions were in the range Ap = 256, . . . , 1024 and At = .01, . . . , .001). 

The CWl is stable in the NLCME context (i.e., stable against transport scale perturbations) and 
according to the dispersive instability criterion (l30l) is also stable for negative values of the disper- 
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Figure 8: CW stability properties for k = 1 and k = 2, with shading indicating instability (dispersive 
instabilities not shown) and dots corresponding to CWl: k = 1, = \, 6 = —'j, and CW2: k = 2, 

p^ = 2,e = l{a = \). 



sion. This is corroborated by the numerical results presented in Fig.9, where the time evolution of 
the spatial norm of the amplitudes and its derivatives, 



is plotted for e = —10"^. The CWl remains stable and the small initial perturbation does not 
grow. Note that the perturbation cannot decay to zero because of the absence of dissipation in the 
system, and thus the value of the spatial derivatives remains finite but small (comparable to the 
size of the initial perturbation, see Fig.9). In contrast, if dispersion is positive, the CWl is unstable 
due to the exponential amplification of the small dispersive scales, see Fig. 10. This growth is not 
a slow time effect; it takes place in the time scale t ~ 1 and does not vanish as £ — > 0, as it can 
be seen in the lower plot in Fig. 10, where the evolution of the spatial derivatives for two smaller 
dispersion values has been added for comparison. The solution that develops is shown in Fig. 1 1 
and it consists of two counter-propagating wavetrains with dispersive wavelength (~ ^/e <^ 1) 
moving at the group velocity; notice how the number of peaks is approximately doubled when the 
dispersion is divided by 4. 

For the CW2 in Fig. 8 the situation is the opposite of that of the CWl. Now the CW2 is 
stable for positive values of the dispersion (the time evolution plot of the norms is omitted in this 
stable case because it is identical to that shown in Fig.9) and the dispersive instability develops, 
producing small dispersive scales all over the domain, for negative dispersion, see Fig. 12. The time 
plot of I l^^l I in Fig. 12 shows the development of the dispersive instability for several dispersion 
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Figure 9: Thick (thin) lines indicate the time evolution of the spatial norm of (A ) and its spatial 
derivative for k = 1 and e = —10^'^. The initial condition is the CWl with a random perturbation of size 
~ 10-3 (a = i). 

0.71 I , , , , , 1 
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0.702 - 

0.7 1 > > > > > 

10 20 30 40 50 , 60 




Figure 10: Thick (thin) lines indicate the time evolution of the spatial norm of A^ (A") and its spatial 
derivative for k = 1 and e = 10~^ (dashed lines correspond to smaller dispersion values). The initial 
condition is the CWl with a random perturbation of size ~ lO"'^ (a = ^). 

values; again, as predicted by the linear stability theory, the dispersive instability exponent goes 
to a nonzero value as the dispersion coefficient is reduced. A space-time representation of the 
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Figure 11: Snapshots of the solution of the NLCMEd at t = 60 (|^+|:thick Hne, |A"|:thin line), with 
parameters as in Fig. 10 and dispersion e = 10"^ (above) and e = 10~^/4 (below). 

solution, once the dispersive scales are well developed, is presented in Fig. 13 where it can be 
clearly appreciated that, in the short time scale t ~ a/RI, the small dispersive scales are just 
transported by the group velocity but they also evolve on a slower time scale t ~ 1 producing a 
very complicated spatio-temporal pattern. 

In conclusion, the main result of this paper is that the standard NLCME formulation fails to 
predict the system dynamics and that, in general, the small material dispersion terms should be 
taken into account, i.e., the NLCMEd eqs. ©-(jT]) should be used, to appropriately describe the 
weakly nonlinear regime of light propagation in a FBG. In particular, the stability of the CW is 
seen to be drastically affected by dispersion. No matter how small is the dispersion or what sign 
it has, there are always stable CW according to the NLCME formulation that are dispersively 
unstable. The destabilization produced by the small dispersive terms is not a higher order, longer 
time effect; it takes place in the time scale of the NLCME and the associated growth rates remain 
of order one as the dispersion coefficient goes to zero. Once the dispersive scales are destabilized, 
they typically spread all over the domain producing a very complicated spatio-temporal dynamics 
(see Fig. 13) that is not captured by the NLCME. This behavior is the result of the competition of 
two effects with different asymptotic order: the dominating advection due to the group velocity and 
the dispersion. This situation generically arises in any extended propagative system with spatial 
reflection symmetry unless some special care is taken to reduce the group velocity. It has been 
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1.5 




Figure 12: Thick (thin) lines indicate the time evolution of the spatial norm of A'^ (A^) and its spatial 
derivative for k = 2 and e = —10"'^ (dashed lines correspond to smaller dispersion values). The initial 
condition is the CW2 with a random perturbation of size ~ 10~^ {a = |). 



A+l \A- 




Figure 13: Space-time representation of the solution in Fig. 12 for two time units after t = 80. 
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previously found in the description of double Hopf bifurcation in dissipative systems [Il7l.ll8h and 



in parametric ally forced surface waves 



There still remains the question of whether one should retain in the envelope equations more 
terms of the form zLae'^A^^^, iPs^A^^xx^ ■ ■ ■ check if they produce new instabilities at even 
higher wavenumbers k ~ kl"^'''^, . . . ^ kl"^^^- This is not necessary as it can be readily 

seen as follows. The two first dominant effects in the evolution of these very short wavelength 
perturbations are transport and dispersion and, proceeding as in section 4, they can be written as 



±ikt+iek'-^t+ikx 



But, because of the effect of the dispersion, there is no possible coupling between with v4^^, 
which is no longer resonant, and hence there is no possible exponential instability (see eqs. (124)) 
and (l25l) ) at these higher wavenumbers. 

As a final remark, it is interesting to mention that there are several recent rigorous proofs 
[a, |20(] that establish that the solutions of the NLCME remain asymptotically close to solutions of 
the original physical equations. These proximity results are not in contradiction with the results of 
this paper because these proofs do not say anything about the stability of the solutions. Therefore, 
when a stable NLCME solution is dispersively unstable, its corresponding close solution of the 
physical equations is also unstable and the system likes to move away from this solution and 
develops dispersive scales, as the NLCMEd correctly predict. 
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